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Abstract 

A new improved study of K~ -proton interactions near threshold is performed using coupled- 
channels dynamics based on the next-to-leading order chiral SU(3) meson-baryon effective La- 
grangian. Accurate constraints are now provided by new high-precision kaonic hydrogen mea- 
surements. Together with threshold branching ratios and scattering data, these constraints per- 
mit an updated analysis of the complex KN and 7rS coupled-channels amplitudes and an im- 
proved determination of the K~p scattering length, including uncertainty estimates. 



Introduction. Within the hierarchy of quark masses in QCD, the strange quark plays 
an intermediate role between "light" and "heavy" . Hadronic systems with strange quarks 
and, in particular, antikaon-nucleon interactions close to threshold are therefore suitable 
testing grounds for investigating the interplay between spontaneous and explicit chiral 
symmetry breaking in low-energy QCD. 

Methods of effective field theory with coupled-channels, based on the chiral SU(3)_r x SU(3) l 
meson-baryon effective Lagrangian, have become a well established framework for deal- 
ing with low-energy KN interactions [1] (see also Ref. [2] for a recent review). However, 
previous applications of such approaches, combining information from earlier kaonic hy- 
drogen measurements [3,4] and older K~p scattering data, were still subject to consid- 
erable uncertainties. The theoretical studies [5,6] gave strong indications for a possible 
inconsistency between the DEAR K~ hydrogen data [4] and the low-energy K~p elastic 
scattering cross section. With the recent appearance of results from the SIDDHARTA 
kaonic hydrogen measurements [7], a new level of accuracy has now been reached that 
permits an improved analysis with updated constraints. The present work describes such 
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an analysis including a new determination of the if~-proton scattering length and im- 
plications for the coupled KN and 7r-hyperon amplitudes below K~p threshold. 

Theoretical framework. The starting point is the chiral SU(3)axSU(3)l meson-baryon 
effective Lagrangian at next-to- leading order (NLO): 

C eff = C M (U)+C$ B {B,U)+C$ B {B,U) , (1) 

where £m(M) with U = u 2 = cxp[i\/2 $//] is the non-linear chiral meson Lagrangian 
incorporating the octet of pseudoscalar Nambu-Goldstone bosons (ir, K, K, rf) in the stan- 
dard 3x3 matrix representation $. At this stage / ~ 86 MeV is the pseudoscalar decay 
constant in the chiral limit. 

The meson-baryon Lagrangian £$ B (B,U) at leading chiral order, £>(p), involves the 
baryon octet (N, A, E, S) collected in the 3x3 matrix field B. The baryon fields couple 
to the mesonic vector and axial vector currents, = (l/2\)\u) d^u + ud^v)] and a p = 
{l/2\)[u^d^u-ud^v)] = -{l/2)u»: 

£^ B =Tr(8(i^-M )B-DB W {a^B}-FB W [a^B^ , (2) 

with the chiral covariant derivative V^B = d^B + i[v^,B]. Here M is the baryon mass in 
the chiral limit. The axial vector coupling constants D and F are determined by neutron 
and hyperon beta decays. At next-to-leading order, 0(p 2 ), the Lagrangian introduces 
several low-energy constants, bi and df 

£<$ B = b Tv(BB) Tr( X+ ) + b D Tt(B{ X+ ,B}) + b F Tr(%+, B]) 

+ di Tr(B {u^ [«", B}}) + d 2 Tr(B [u„, [«", B}}) 

+ d 3 Tv(Bu tl )Tr(Bu^ + d 4 Tr(BB) Tr( U/i u") , (3) 

where x+ = 2 Bo (uA4u + Mu^ is the symmetry breaking term with Bo representing 
the magnitude of the chiral condensate divided by / 2 , and M. — diag(m u , m d , m s ) is the 
quark mass matrix. At tree level in chiral perturbation theory (ChPT), the constants 
bo, bo and b F are constrained by the baryon octet masses (their splittings and shifts 
from the chiral- limit baryon mass M ). Note however that the present analysis goes well 
beyond tree level so that these constants need not be identical to the ones from ChPT. 
They are renormalized by loop effects taken to all orders. 

We recall that C^b °f ^q. (2) generates the leading Tomozawa- Weinberg (TW) terms 
of the interactions between the meson and baryon octets (Fig. 1(a)). It also generates 
meson-baryon direct and crossed Born terms (Fig. 1(b) and (c)) upon iteration of the 
pseudovector derivative coupling vertices proportional to D and F. The NLO terms 

(2) 

(Fig. 1(d)) derived from £ MB of Eq. (3) involve, apart from the bi coefficients, four low- 
energy constants dj that will be varied freely to achieve a best fit to the available KN 
threshold and scattering data. 

The driving meson-baryon interactions of Fig.l, derived from the NLO effective La- 
grangian (2) and (3), serve as input interaction kernel (denoted by Vy) for the coupled- 
channels Bethe-Salpeter equations connecting meson-baryon channels i and j. We use 
the full set of ten strangeness S = — 1 channels with index assignments i — 1, . . . , 10, 
provided by the baryon and pseudoscalar meson octets and numbered in this order: 
i = K~p, K°n, 7T°A, 7r°S°, Tr+ST, tt-£+, t]A, ??£ , K+E~ , K°~°. The % depend on 
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(b) 



(c) 



(d) 



Fig. 1. Driving interactions generating the mcson-baryon coupled-channels amplitudes: (a) leading order 
(Tomozawa- Weinberg) term; (b) and (c): direct and crossed Born terms; (d) NLO terms. Dashed lines 
represent SU(3) pseudoscalar octet mesons, solid lines refer to members of the baryon octet. 

the meson-baryon center-of-mass energy, i/s, the scattering angles, f2 = {9, if}, and the 
baryon spin degrees of freedom, <7j ; j. We concentrate on s-wave driving terms, Vij(y/s) = 
(l/87r) / dfi ^(-^/s, fi, cr), summarized and explicitly listed in Refs. [2,5,6]. For ex- 
ample, the leading order Tomozawa- Weinberg term is simply 1 



V; 



(TW) 



(V5) = -^(2V5-M i 



Mrfy/iMi + EiXMi + Ei) 



(4) 



where Mj and = M 2 + q 2 are the baryon mass and energy in channel i, with 
(\i the center-of-mass momentum in that channel. The constants C'ij are determined by 
SU(3) Clebsch-Gordan coefficients and given in Refs. [1,5,6]. The s-wave coupled-channels 
T— matrix with elements is found by solving the matrix integral equations 



V + V G T = (V- 



G)" 



(5) 



Here G is the Green function matrix. Its elements Gij = Gi(q 2 ) Sij are the meson-baryon 
loop functions, 



G l (q 2 ) = 



d 4 k 



(2tt) 4 [(q -k) 2 - Mf + ie] (fc 2 - mf + ie) ' 
evaluated in each channel i using dimensional regularization: 



(6) 



Gi(y/s) = ai(fj,) + 
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(7) 



where mi is the meson mass in channel i. The subtraction constants ai(fi) act as renor- 
malization parameters at a scale /x. They cancel the scale dependent chiral logarithms 
and make sure that the calculated observables are scale invariant. 

Observables. Forward scattering amplitudes and cross sections are given as 



f 



Stt^ Tij 



and a lJ (y/s) 



16-7TS 



The K p scattering length is a(K p) = lim t /, f(K p — > K p) = 
Further observables of interest are the threshold branching ratios 



(8) 

h 1 {sfs = m K - +M P ). 



1 The normalization convention used here is the same as in Ref. [5], with dimensionless Vij. It differs 
from the one used in Ref. [2] by a factor MiMj. 



3 



r(tf-f>^7T+£-) 0-51 D r(X-p^^ A) CT31 
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r(if p — > 7r E+) (76i ' ™ r(if p — > neutral states) (T31 + CT41 



(9) 



r(X p — > all inelastic channels) 031 + 0-41 + 051 + ogi 



with all partial cross sections (Jij taken at K~p threshold. The (ancient but accurate) 
empirical values of these branching ratios [8] are listed in Table 1. The energy shift and 
width of the Is state of kaonic hydrogen are related to the complex K~p scattering length 
as 



with the K p reduced mass, fi r = muMp/ (mx + M p ), and including important second 
order corrections [9]. We use the accurate SIDDHARTA measurements [7]: 

AE = 283 ± 36(stat) ± 6(syst) eV , r = 541 ± 89(stat) ± 22(syst) eV . 

The available data base is completed by the collection of (less accurate) scattering cross 
sections [10] (see Fig. 2). We do not include measured irT, mass spectra in the fitting 
procedure itself but rather generate them as "predictions" from our coupled-channels 
calculations. 

Results and discussion. Using the unitary coupled-channels method just described, the 
basic aim of the present work is to establish a much improved input set for chiral SU(3) 
dynamics, by systematic comparison with a variety of empirical data and with special 
focus on the new constraints provided by the recent kaonic hydrogen measurements [7]. 
A detailed uncertainty analysis is performed. It will be demonstrated that previous un- 
certainty measures [5,6] can be reduced considerably. 

We have carried out x 2 fits to the empirical data set in several consecutive steps: first 
starting with the leading order (TW) terms, then adding direct and crossed Born terms, 
and finally using the complete NLO effective Lagrangian. The results are summarized in 
Table 1. All calculations have been performed using empirical meson and baryon masses. 
This implies in particular that those parts of the NLO parameters bo, bp and bp respon- 
sible for shifting the baryon octet masses from their chiral limit, Mq, to their physical 
values, are already taken care of. The remaining renormalized parameters, denoted by 
bo, br> and bp, are then expected to be considerably smaller in magnitude than the ones 
usually quoted in tree-level chiral perturbation theory. Similar renormalization arguments 
imply that the pseudoscalar meson decay constants should be chosen at or close to their 
physical values [11], 



It turns out that best fit results can indeed be achieved with these physical decay con- 
stants as inputs. This is a non-trivial observation, as previous calculations have com- 
monly used an average decay constant as a mere fit parameter, irrespective of physical 
constraints. 

With the TW terms alone a reasonable overall fit (with x 2 /d.o.f. = 1.12) can be reached 
but the kaonic hydrogen energy shift comes out too large (AE = 373 eV) and some of 
the subtraction constants dj in Eq. (7), especially those in the irA and ijT, channels, 
exceed their expected "natural" values ~ 10~ 2 by more than an order of magnitude [12]. 



AE - if/2 = -2a 3 ^ 2 a(K~p) [l + 2a fi r (1 - In a) a(K~p)] 



(10) 



/„ = 92.4 MeV , f K = (1.19 ± 0.01) /, 



/, = (1.30 ±0.05) f n ■ 



(11) 
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TW 




TWB 


NLO 


Experiment 


AE [cV] 


373 




377 


306 


283 ± 36 ± 6 


[7] 


r [ev] 


495 




514 


591 


541 ± 89 ± 22 


[7] 
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2.36 




2.36 


2.37 


2.36 ±0.04 


[8] 


Rn 


0.20 




0.19 


0.19 


0.189 ±0.015 


[8] 


Rc 


0.66 




0.66 


0.66 


0.664 ± 0.011 


[8] 


X 2 /d.o.f 


1.12 




1.15 


0.96 




pole positions 


1422 - 


16i 


1421 - 17i 


1424 - 26i 






[McV] 


1384- 


90i 


1385 - 105i 


1381 - 81i 







Table 1 

Results of the systematic x 2 analysis using leading order (TW) plus Born terms (TWB) and full NLO 
schemes. Shown are the energy shift and width of the Is state of the kaonic hydrogen (AE and T), 
threshold branching ratios (7, R„ and R c ), x 2 /d.o.f of the fit, and the pole positions of the isospin 7 = 
amplitude in the KN-nT, region. 

This clearly indicates the necessity of including higher order terms in the interaction 
kernel Vjj . It also emphasizes the important role of the accurate kaonic hydrogen data in 
providing sensitive constraints. 

The additional inclusion of direct and crossed meson-baryon Born terms does not 
change AE and x 2 /d.o.f. in any significant way. It nonetheless improves the situation 
considerably since the subtraction constants a, now come down to their expected "nat- 
ural" sizes. 

The best fit (with x 2 /d.o.f. = 0.96) is achieved when incorporating NLO terms in the 
calculations. The inputs used are: the decay constants f n = 92.4 MeV, Jk = 110.0 MeV, 
f v = 118.8 MeV, and axial vector couplings D = 0.80, F = 0.46 (i.e. g A = D + F = 1.26); 
subtraction constants at a renormalization scale fj, = 1 GeV (all in units of 10 -3 ): a\ = 
a 2 = —2.38, 03 = —16.57, 04 = a 5 = a 6 = 4.35, aj = —0.01, as = 1.90, ag = a w = 
15.83; and NLO parameters (in units of lO" 1 GeV" 1 ): b = -0.48, b D = 0.05, b F = 
0.40, di = 0.86, d 2 = -1.06, d 3 = 0.92, d 4 = 0.64. Within the set of altogether 
"natural" -sized constants a* the relative importance of the KE channels involving double- 
strangeness exchange is worth mentioning. 

As seen in Table 1, the results are in excellent agreement with threshold data. The 
same input reproduces the whole set of K~p cross section measurements as shown in 
Fig. 2 (Coulomb interaction effects are included in the diagonal K~p — > K~p channel as 
in Rcf. [5]). A systematic uncertainty analysis has been performed by varying the input 
parameters within the range permitted by the uncertainty measures of the experimen- 
tal data. A detailed description of this analysis will be given in a longer forthcoming 
paper [13]. 

Equipped with the best fit to the observables at K~p threshold and above, an opti- 
mized prediction for the subthreshold extrapolation of the complex s-wavc K~p — > K~p 
amplitude can now be given. The result is shown in Fig. 3, including again a conservative 
uncertainty estimate. The real and imaginary parts of this amplitude display as expected 
the A(1405) resonance as a quasibound KN (I — 0) state embedded in the irT, contin- 
uum. The present NLO calculation confirms the two-poles scenario [14,15] of the coupled 
K~p o 7rE system. Using the best-fit input, the resulting locations of the two poles 
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Fig. 2. Calculated K~p elastic, charge exchange and strangeness exchange cross sections as function of 
K~ laboratory momentum, compared with experimental data [10]. The solid curves represent best fits 
of the full NLO calculations to the complete data base including threshold observables. Shaded areas 
give an impression of uncertainties. 

in the complex energy plane are as follows: "upper" pole ( A'.V-dominated) : 1424 - i 26 
MeV; "lower" pole (7rS-dominated) : 1381 - i 81 MeV. Unlike previously found patterns 
in which the location of the lower pole has been subject to large model uncertainties, 
the pole positions now remain remarkably stable with respect to changes of the input. 
The shift of the real parts of both these pole positions from the "TW" and "TW + Born 
terms" steps to the full NLO calculation is less than 5 MeV. The corresponding change 
in the imaginary parts is only slightly larger (between about 10 and 20 MeV). 

The K~p scattering length, a(K~p), deduced from the kaonic hydrogen measure- 
ments [7] and with inclusion of Coulomb corrections (see Eq.(10)) is: 

Re a(K'p) = -0.65 ±0.10 fm , Im a(K~p) = 0.81 ± 0.15 fin , (12) 

with an error estimate based on the uncertainties assigned to the measured kaonic hy- 
drogen energy shift and width. Our best fit NLO result, a(K~p) = —0.70 + i 0.89 fm, 
is perfectly consistent with Eq. (12). Note that this new determination of the K~p scat- 
tering length has shifted quite significantly in the value of Re a(K~p) from previous 
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V» [MeV] [ MeV ] 

Fig. 3. Real part (left) and imaginary part (right) of the K~p — > K~p forward scattering amplitude 
extrapolated to the subthreshold region. The empirical real and imaginary parts of the K~ p scattering 
length deduced from the recent kaonic hydrogen measurement (SIDDHARTA [7]) are indicated by the 
dots including statistical and systematic errors. The shaded uncertainty bands are explained in the text. 

ones [5,6,16,17], mainly because of the new constraints from the much improved SID- 
DHARTA data, 

Summary. Given the significantly more accurate constraints from the new kaonic hy- 
drogen measurements, an improved theory of low-energy antikaon-nucleon interactions 
on the basis of chiral SU(3) effective field theory with coupled-channels is now at hand. 
The results and conclusions are summarized as follows: 

1) Kaonic hydrogen data are now consistent with low-energy K~p elastic, charge ex- 
change and strangeness exchange cross sections. 

2) The present next-to-leading-order (NLO) analysis indicates a well convergent hierar- 
chy of driving terms for coupled-channels dynamics derived from the chiral SU(3) meson- 
baryon effective Lagrangian. The Tomozawa- Weinberg terms dominate, Born terms are 
significant, while a "best fit" nonetheless requires NLO contributions, though with rela- 
tively small coefficients. 

3) A new, more accurate determination of the K~p scattering length has been pre- 
sented. 

4) As an important result of the present analysis, the best NLO fit systematically 
prefers physical values of the decay constants f v , fx and f v . This is actually a (successful) 
test of consistency, given that physical meson and baryon masses and empirical baryon 
axial vector coupling constants are used as input. 

5) The two-poles scenario of KN and 7rE coupled-channels dynamics is re-confirmed. 
The predictions for the pole positions in the complex energy plane have been sharpened. 

6) Uncertainties in the subthreshold extrapolation of the KN amplitude arc reduced 
as compared to previous work. 

A more detailed presentation including further results and predictions, e.g. on 7r£ in- 
variant mass spectra and the K~n scattering length, is in preparation. 
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